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ABSTRACT 


Performance improvement of a carbon dioxide gas-dynamic 
laser, due to regeneration, was studied theoretically. The 
study determined to what extent preheating the nitrogen, 
before mixing it with the carbon dixoide and water, would 
effect available laser output energy. Combustion chamber 
pressure and nitrogen inlet temperature were incremented 
by 100 psia and 300 K respectively. Due to the fact that 
little dissociation occurred at temperatures studied, no 
pressure dependence was found for the throat or combustion 
chamber temperatures. A method of preheating or regenerating 
the nitrogen is to use a gas-to-gas heat exchanger at the 
exhaust of the diffuser. Small signal gain was improved 
maem 0.35 to 0.78 mt. Limiting power extraction was 
increased by 42%. Optimum nitrogen temperatures were 
determined to be 1350 K for small signal gain and 2400 K 
for optimum power extraction. For typical operating con- 
ditions the fractional pressure loss across the heat 


exchanger was found to be 0.1 atmospheres. 
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exchanger total heat transfer area on one side. (£t7) 
exchanger free flow area. (£t7) 
Avagadro's Number (6.02 x Mae Mole“) 


constant pressure specific heat 
Eanmnang £riction factor 
exchanger flow stream mass velocity [lbs/ (Hr-Ft*) 3 


gravitational constant [32.2 (£t/sec*) (>) ] 
m 


Piranck S Constant [6.62 "x ior? J-Sec] 
saturation intensity (Kw/m*) 
output beam intensity (Kw /m-) 
rotational quantum number, 
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Boltzimann's constant [1.38 x 10 2 J/K] 


Gone EaActIOn LOSS COGCELrIclent for flow at heat 
exchanger entrance. 


expansion loss coefficient for flow at heat 
exchanger exit. 


Mach number 

mass of one molecule of xx. 

mass flow rate of xx. 

pressure in cavity (region 4 in Figure 1) psia 


rate of heat transfer from hot gases to cold 
nitrogen. 


universal gas constant 
nitrogen temperature (K) 


combustion chamber temperature (K) 
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mixing chamber temperature T, = Ty (K) 
nozzle throat temperature (K) 
translational-rotational temperature in the cavity (K) 


upstream specific volume of gases with respect to 
heat exchanger. 


downstream specific volume of gases with respect to 
heat exchanger. 
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Z 
small signal gain coefficient (m 


mean specific volume = 
=) 
C/C,, = 1.4 for this mixture of gases 
effective optical cross section 


ratio of free flow area to frontal area 


optical broadening collisional cross-section of 
Sa (om ) 


characteristic vibration temperature of CO, for 


vi» (1995 K) é 
characteristic vibration temperature of co. hor 
Vo. (960 K) 


characteristic vibration temperature of No (S300) 


characteristic rotational temperature of CO 
(2.88 K) 


motemirtaction Of xx. 

density of gases 

laser frequency (2.83 X oe) Hz 
optical collision frequency 
radiative lifetime (5.38 sec) 
MeGucea Mass Of Xx 
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wavelength of laser radiation (10.6 X 10 mM) 
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cen tRODUCT ION 


The carbon dioxide gas dynamic laser has been considered 
a prime candidate for many high energy laser programs. 

One problem that has always plagued the high energy laser 
program is that of excessive weight and volume. These 
presently are the factors that are impeding the development 
of a tactical Airborne laser weapon system. A heat exchanger 
can be placed at the exhaust of the diffuser to preheat the 
nitrogen gas. This system significantly reduces the weight 
and volume of a given gas dynamic system while maintaining 

a constant level of energy output. 

The combustion chamber properties and nozzle throat 
conditions were determined by using a computer program 
developed by Sanford Gordon and Bonnie McBride at the NASA 
Lewis Research Center in Cleveland, Ohio. The Program and 
a detailed description of its use are contained in reference 
{1}. Although the program was written with a rocket engine 
analysis in mind, it is adequate for computing equilibrium 
and frozen composition conditions in combustion gases used 
in gas dynamic lasers. 

After the computer determined conditions up to and 
including the nozzle throat, hand calculations were made 
using methods specified in reference [2]. According to 


calculations, the maximum available energy can be increased 





by a factor of 5 to 10. The small signal gain as a function 
of nitrogen temperature reached a maximum value of 0.773. 
This maximum indicates a best operating temperature. 

The diffuser is used tc recover the pressure in the flow 
stream and must be matched to the heat exchanger. 

A gas to gas heat exchanger is used to insure that the 
proper temperature difference is attained from the cold 
nitrogen side to the hot nitrogen side (see Figure 1). The 
pressure must not be allowed to decrease to atmospheric 
pressure at mean sea level. This can be done by use of 


procedures outlined in references [3], [4], and [5]. 
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Pas ORY 


fee BASIC THEORY OF A GAS DYNAMIC LASER AMPLIFIER 
Carbon monoxide, hydrogen and oxygen are fed into the 
combustion chamber and ignited. The amount of each reactant 


is determined by the following reaction; 


te COmwenem toons JCO- + HAO + YON 


JON, 2 2 2 2 2 


After the reactant products reach equilibrium, nitrogen is 
mixed with the carbon dioxide and water in the mixing chamber. 
This yields a molar mixture of 90% nitrogen, 1% water and 
9% carbon dioxide. This is the composition required for 
Operation of a gas dynamic laser amplifier. It is assumed 
that the ratio C./ oy Sempenecemixtunue 1c 1.4 while the 
Vibrational mode of carbon dioxide is frozen out at T, > 

The nozzle array is designed to attain Mach 1.0 at the 
throat (region 3 in Figure 1) and Mach 4.0 or greater inside 
the cavity (region 4 in Figure 1). For purposes of this 
thesis, all calculations have been done assuming that a 
Mach number of 4.0 exists inside the amplifier cavity. The 
technique used to determine the area ratios is described in 
detail in reference [6]. Since these calculations are well 


documented, only some of the final results will be mentioned. 


dg 





From the equilibrium combustion temperature, the throat and 
cavity temperatures can be determined and the small signal 
gain can be calculated. 

The small signal gain is applicable only inside the 
cavity (region 4 in Figure 1). Assuming that the oscillator 
and cavity are matched at a wavelength of 10.6 micrometers, 


the output beam intensity is given by Equation l. 
ae ae Q ot (1) 


a is the gain coefficient and 2 is the overall beam pathlength 


inside the cavity. 2& is found from Equation 2, 
£ = NL (2) 


N is the number of times that the beam traverses the width 
L of the cavity (see Figure 1). 

A measure of the power that can be extracted, defined 
in Equation 5.23 of reference [2] as the limiting power 


extraction, is aot. where I, is defined by Equation 3. 


eC — 5) 
a) oe Le (3) 
re) Ey 
'e, TItT A) E32, 


The maximum available energy in the excited state of the 


gases is deScribed in reference [2] and can be calculated 


iL 





with Equation 4. Table II shows the result of calculations 


at incremental temperatures. 





ae il hv 
Emax ~ = 3380] oe 
eee Te -l 


m is defined as the average mass of the gas molecules. m 


is calculated with Equation 5, 


y M io) some! 
7" SO. ED, HO HO N,N, 7 


A 
V 


y is the mole fraction of the subscripted quantity. 

As a basis for determining the change in limiting power 
output, small signal gain and maximum available energy, an 
assumption was made that all gases were at a room temperature 
of 300 K, prior to combustion. The nitrogen temperature was 
then incremented and new calculations made at each value of 
nitrogen temperature. These new calculations are compared 


to those at 300 K in tabular and graphical form. 


B, SMALL SIGNAL GAIN 

Reference [2] outlines procedures for determining the 
Small signal gain coefficient. The overall gain I/t, may 
be determined from Equation 1 after a has been found from 
Equation 6. Pages 5-12 of reference [2] shows the relation 


between c and ao. 
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~J(J+1) 0, 


2 
a = 2 Oo (2a+1) 2 “4 n(001)-n(100), yy 
fe) 47 ie T7257. Yao N CO. 


(6) 
The population inversion is determined by the quantity in 
brackets in Equation 6. Population inversion is related to 
cavity and throat temperature by Equation 7. It is assumed 
that the vibrational temperatures of nitrogen and carbon 
dioxide are frozen at T3- Appendix B is a derivation for 


oy in terms of temperature and number density. 


-6_/T =i I 
po(001)-n (100) | _ -e i e 4 
aie je (Se 5 (7) 
CO. VIB 


Here, the temperature T. is frozen at the nozzle throat and 
the system is operating so that the threshold condition is 
exceeded. That is, 6/8, < T/T). This condition insures 
that the small signal gain will be positive. Q.., is related 


vib 
to throat and cavity temperatures by Equation 8. 


-6,/T -6,/T -6,/T 
6) ee oles 1 Ay mt ce 2 25, 7G Files 5 3,-1 


(8) 
C. MAXIMUM AVAILABLE ENERGY 

The maximum available energy is dependent upon the 
wavelength of the laser and the vibrational temperature of 


Eiencwansitren level for stimilated emission. oe can be 


Ae, 





calculated using Equation 4 or Equation 9. Equation 9 is 
the same aS Equation 4 except that 5 has been substituted 
mer Vv. 


6,/T, oi 


he 5 a 
max (e = 1) (9) 


ry 
I 
Sy ft 


D, HEAT EXCHANGER 

The heat exchanger is the key to success of this system. 
It takes energy from the exhaust gases at the diffuser and 
transfers it to the nitrogen gas. This process allows the 
nitrogen to be heated to an optimum temperature. Equation 
13 deScribes the pressure drop across the exchanger. Since 
this system is to operate at mean sea level and higher, the 
minimum pressure on the downstream side (region 7 of Figure 
1) of the exchanger must be greater than 1 atmosphere. 

Assuming that the nitrogen enters the exchanger at 300 K 


and exits the exchanger at a temperature of T,, (optimal), 


N 


the total heat transfered is defined in reference [3] as 


Eoeuation 10. 
aN = mM CoN [T., (HOT - Ty (COLD) J (10) 


Se must equal the heat given up by the mixture flowing thrcugh 


the heat exchanger. This is described by Equation ll. 


Gy = q =m Cc [T,(HO) - T_ (COLD) ] (11) 
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Therefore, Equations 10 and 11 can be combined to produce 
Equation 12 which describes the temperature drop across the 


heat exchanger. 


_ Ny, (pn 
Ts — Ts = ae ae [T., (HOT - Ty (COLD) ] (152)) 
m pm 


Here, Te. is defined as the upstream temperature, Ta. (HOt), 


and T. is the downstream temperature, Th. (elec yee 


5 
A more detailed analysis is in Appendix A. Equation 13 


is used in reference [4] to deScribe Ap. 


2 
Gvy 


= 2 


V V V 
*) 42 (22-1) +f OM - (1-0°-K,) 41 
1. eG i 
(13) 


Ly 





Pile CALCULATIONS 


A. THE GAS DYNAMIC LASER CAVITY 
Calculations for the combustion chamber and nozzle 
throat were computed by using the NASA Lewis Computer 
Program, reference [lJ. The remaining calculations were 
done by hand uSing the isentropic tables in reference [6]. 
For conditions in region 4, shown in Figure 1, Appendix 
A of reference [6] was used. At M = 4.0 the following 


relations are found: 


ele f 
Pp y-1 
== (1 + yl M*) = 0.6586 x 10 Z 
2 
i 
T —J1 
ee, ve oe = 
rT (1 + saat ) = 0.2381 
A 
A 4 
See = 10,72 
— ; 
Arhroat 


Conditions in region 4 (Figure 1) can be determined from 
the following relationships: 
2 


Pe) Oe o> S566% LO. ) 2 (14) 


(OR zoc1) Ty (15) 


FF 
Il 
1: 
=. 
il 
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AY, 


— = (10.72) A 
throat 


By = 1 ) (16) 


* hroat throat 


One model of the starting process uses a normal shock 
which develops inside the cavity (region 4 of Figure l). 
Therefore, the diffuser must be able to swallow this shock 
before the system will be stable enough to operate effec- 
tively. From Appendix B of reference [6] at M, = 4.0, the 
Mach number of the flow behind the normal shock is 0.4350. 
Before the diffuser can swallow the shock, the flow must be 
accelerated to Mach 1 at the diffuser throat. Therefore 


from Appendix A of reference[6] at M = .4350, the following 


ratio is found: 


A 
eee = 4S (17) 
* 
AX Phe eaeer 


Diffuser conditions change after the shock has been swailowed. 
Therefore, the only noteworthy condition is the fact that 
the ratio of diffuser throat area to nozzle throat area can 


now be determined. 


* 
Agi ffuser a A/S 2 JUG ese 7°45 (18) 
Athroat A/AS 1. 


By the methods described in Section 3.4 of reference [6], 


the Mach number at the diffuser throat is now 3.55 and the 


Eg 





normal shock is positioned just past the diffuser throat 


near region 6 of Figure l. 


B. SMALL SIGNAL GAIN COEFFICIENT 
The small signal gain coefficient must be determined 
before the output intensity can be calculated. Combining 
Equations 6, 7, and 8 yields Equation 19, Equation 19 is 
explicit in terms of known throat and cavity temperatures. 
The constant 1.6 x Lye is the result of inserting 
6 


2 = - = = ; -14 
a oero.6 x 10 m, J = 20, oa 2h este) Ne ve Seige ee ILO) N VT," 


and +t = 5.38 sec. 





-239.4 -3380  -1995 
2 = v4 e = - e “4 
Od, = [1.6%10°) [———y79) (=p 995-9 0 rr 
(T4) ir, 7 = 
Gee 9 Gigs “Galea ee 
(19) 


Table 1 gives the values of nitrogen, throat and cavity 
temperatures used in calculating a, at the specified nitrogen 
temperatures. Figure 2 is a graph of a, versus temperature 
of the nitrogen. It is shown that an operating nitrogen 
temperature of approximately 1300 K will maximize the small 
Signal gain at 0.773. See Table II for values of a, at 


specified nitrogen temperatures. 
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C. MAXIMUM AVAILABLE ENERGY 
The maximum available energy can be easily calculated 
by using Equation 4 after m has been determined by using 


Equation 5, 





mee) (44) + (,07) (18) + (.90) (28) _ -23 gm 
a D ee Me oracuic 
6.02 x 10 
3380 
T 
a 5 3 -l Joules 
Eee (3-283) x 10) (e iL Ten 


These results are listed in Table II and Figure (3) is a 


prot of eee versus nitrogen temperature. 


D. LIMITED POWER EXTRACTION 

The saturation beam intensity, as defined by Equation 3, 
is a valuable aid in approximating the limiting power extrac- 
tion. Equation 3 must be used in conjunction with the 


following equations: 


-~6._/T 6./T -1 8 ST ALR -i -l 
_ 4 Ly gee 2 DoeA 
e =e [(e -1) + 2 o, (e ~ 1) ] (20) 
~6._/T 

ee = (a, (21) 
; -6 J (J+1) 
r > a 

rp = + 98  (23tl) , 4 (22) 


o064T ie 170m 


ON 





Senpining EGuations 3, 8, 20, 21 and 22 yields Equation 23. 


a ae 
T i TT 
po 8378) re 4)F re 4) 77 (1-e 3) 7+ ve 
oO E 2 
(1 +14) iL ho 204+1 Ue 
3 aT TVn T7265 


The product of Equation 19 and Equation 23 is Equation 24. 





3 a 
fy Aly 
a I §6.-6 3 4 
mee = [10Y.5 | ie) {——_— S—_____ (24) 
A 2 : Ee _- Pattie) 
3 E. 14 4 4 


TP, and TAP, can be found in 


Figures 2.8 and 2.10 of reference [2]. Table III shows the 


oe ; 
Ot a/ Py 1s in Mks units. 


calculated results from Equation 24 and Figure 4 is a pilot 


opi a 1 /P versus nitrogen temperature. 


BE. HEAT EXCHANGER 

Using the optimum temperature as shown in Figure 4 
the hot nitrogen temperature is specified as 2400 K. This 
gives a temperature difference in the nitrogen of 2100 K. 
An energy balance will produce an equation relating the 
temperature difference from the hot to the cold side of the 


heat exchanger. From equation 10 and C equal to 06.39, 


pN 
the rate of heat transfer is 


22 
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e a e = e BTU 
q = my CON (AT) m, (819) sp——ae=. 
Reference [7] contains values for oo From Equation 12, 
the temperature difference is 


Cc N 
T - T. = (1803) oa 2 1617 K 


> 7 
pm 
A more detailed analysis of the heat exchanger is in 


Appendix [A]. 


F., APPROXIMATE WEIGHT AND VOLUME 

Assuming One-Second bursts of energy in the output, the 
amount of fuel needed for any pre-determined number of 
bursts can be calculated eaSiest at the nozzle throat where 
gas properties are well defined. From reference [6], the 


mass flow rate is 


re 
il 


VAP 


where U is the mass velocity, A is the area of the nozzle 


and p is the density of the gases, 


The limiting power per unit mass is 


ZS 





where m = (1 sec) ma If W increases by some factor A and 
m remains constant then aot. must have increased by A also. 
It is desired to maintain W constant and adjust the mass 


necessary to compensate for this increase in output. 


(1 + A) a1, 








W(1l + A) = 
m 
ae Se m 
= =) m' = 
fi i (145) 
(i+A) 


This means that the mass neceSsary to maintain a constant 

3 A Pe 
output can be decreased by a factor of TitAy ° Calculations 
have shown that A = 0.42 and the fuel necessary is only 703% 


of the fuel used without a regenerator. Therefore a savings 


in fuel weight and volume of 30% has been found. 
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IV. RECOMMENDATIONS AND CONCLUSIONS 


It has been shown that pre-heating the nitrogen gas 
prior to mixing it with the combusted products of hydrogen, 
oxygen, and carbon monoxide increases the gain of a carbon 
dioxide gas dynamic laser from 0.35 he te OmOn 07 mo, The 
limiting power extracted reaches a maximum at a temperature 
of 2400 K. This produces an optimum operating temperature 
and the limiting power extracted can be increased by a factor 
of 1.42. This implies an increase of 42% in output power 
while maintaining a constant mass flow rate. This means that 
a constant level of output power can be maintained while 
reducing the amount of fuels needed by 30%, This decreases 
the amount of storage volume by 30%. 

It is recommended that a further study be conducted to 
determine the effects that a heat exchanger would have on 
other fuels to produce a 10.6 micrometer wavelength with a 


possible extension into other wavelengths. 
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TEMPERATURES AND PRESSURES OF SYSTEM 
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TABLE I (Continued) 


AGH 


) 
500_ macs 3 
600 Sy sied 
700 |_ 4.610 
1450 1975 1730 470 wecG0 ih esos. 
200 
| 1000 


O10 
600 
1500 2014 480 inn 
1800 2246 1991 : 
2100 2463 2204 
q _ | 
2700 2845 2590 677 aati 510 
F900 | 5,927 
1000 6 586 
500 
600 
3000 700 4.610 
3017 2760 718 as 610 
900 5,927 
1000 6, 586 


20 
















TABLE II 
COMPUTED VALUES FOR aay & om 


max 


K K/KG) «CE + KJ/1bm a, (m=") 
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TABLE III 
COMPUTED VALUES OF ap1,/P4 oe 
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APPENDIX A 


THEORY OF GAS-TO-GAS HEAT EXCHANGER 


As stated on page 21 of reference [4], the terms 
involving exchanger input and output flow are negligible. 


An inspection of Equation 13, shows that the assumption 


< 


V V 
m Z Z 2 Z 


iL ae 


a 


f a 
S 
can be used to develop Equation 4', K t1-0" 1s associated 
Vv 
with the upstream flow while (1-6*-K_).* is associated with 
iI 


the downstream flow. 


2 


Gov V Vv 
a5 it Zn A om 


By definition: 


—_— 1 2 
m y e 
Therefore: 
2 
Gvv Vv 
— Ih 2 ears fA 
cme ry erm e  *)! 
& Cc 


Now, the heat transfer area is needed. The derivation 
is based on the counterflow heat exchanger, from reference 


(7), shown below. 
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Ty (HOT) 


© @ 


———1 > Se 


Ty (COLD) 


The temperature variation as a function of distance 


through the heat exchanger is shown below, 





is 
AT 
T (RK) 
T, (HOT) T. 
Ty (COLD) 
DISTANCE 

Now: 

. BTU 

q = hAAT (sao) 

AT = logarithmic mean temperature difference. 
h = heat transfer coefficient, 
A = transfer area 


Before A can be determined, AT is needed. Therefore: 


a(AT) = dT +. dT,, 





where: 


Now: 


Therefore: 


Also: 


Therefore: 





qv = change in hot gas temperature 
qty = change in cold gas temperature 
Sigi = 2G el GUE eee cue 
a(AT) = oe oc ee 
“Cn Mn Cy My 
dq = hAT @A 
cea = h(—~ - —+ +) aa 
e Gu i Cc. om 
n "N m om 


Integration of this equation leads to 


To-T,, (COLD) 


1g ler THOT! 
Te ia OT 


il 


Cm 
mm 


if 


Os Ms 








hA ( ) 
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As shown previously 


hA = q/AT 


and 
: M Cm (Ts5-T7) — MyCy (Ty (HOT) - Ty (COLD) ] 
q/VI = 
AT a AT 
Now: 
T5-T, (COLD) _ T. (HOT) ~ Ty (COLD) _ T. + Ts 
Cae iory saa ~ ° Ar 
5 N m 
and 
_ (T.. (HOT) = T,, (COLD) - qT. + To] 
AT = . 
T - T (COLD) 
i N 
ln a 
5 N (HOT) 
Now: 
a m_C,,(Te-T>) _ Ps U, A.C (fT. - To) 


hAtT hat 
m m 


Amd the ratio A/A, can be found as 


[ poUeC (T.-TW) ] / (hav. | 


Si 





Assuming a 10% pressure drop across the exchanger, Pe. 
equal to 1.1 atmospheres and T, = 4500 °R, the Reynolds 
number can be calculated. 
DU 
O 


R = 
e 


5 Ps 
u 


Assuming a tube diameter of .1 ft, Uz = 988 (ft/sec), 


: 4 STE _ =) Seg 
o> B10 2 x 10 wre) and Le = 1,09 x 10 ft-sec’ 


Re = De. x 10° 


From Figure 6-9 of reference [3] 


oa Z 
a Sf oe oO 


k = K ef LAr 


ref 


where k and T are from Table a-5 of reference [3]. 
ia ie ref 


4500 _ 05 BDU 


h = 336 BTU/(hr £t* °R) 
So, °C 
— == (tf, - fT) = 12 
Cc hat. 
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Vv 
With these quantities and —~ equal EOuleonand £ = 0.02 
1 


‘(from Figure 6-3 of reference [4] and c*v, =o. 36) x 104 


1b*/(lbm fete SG67)) ¢ 


2 
sie 1b 
AP x 3g = 0.1 atmospheres = 212 > 





REPENDEX S 


DERIVATION OF es 


From Equation 3.7 of reference (2], Vo is 


oa C.. ite Oa, CO Cem Wie. 0... Ci] 
a a CO. CO. CO. N. No N. HO HO HO 


wy = the mole fraction 
0 is the optical broadening collisional cross 
section. 


C, is the relative mean speed. 


From pages 3-5 of reference [2] 


fe) = 1.3 x 10714 erie 
co. 
ee on] ven: 
N 

2 
Oo = ,38 x 10774 oRe 
HO 
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From page 3.4 of reference [2] 


= _ (tS eas 


on 4 
al = av 
where 
y= RCO, for CO. Gol sions 
elo Me. 
y= ———-_ for HO and N. collisions. 
M +M 
CO me 
2 
Now: 
= [8 K Ty M 
comme 7 M age V4 Sec 
2 CO. 
= 8 RT,’ 7 
e = Ee 44,2 T,' — 
N 4 sec 
2 TT Hy 
2 
_ [? R T," M 
CH one = 40 Vt4 sec 
2 | T. 
H,O 
2 
Therefore: 


v= N fT x 107741 (21) (.09) (1.3) + (44.2) (.9) (.87) 


PO) GeO) (5.38). | 


3.95 x 10-° N (2, 


Al 
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